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Abstract 



In this article we derive Fokker - Planck equation for incompressible fluid and investigate its prop- 



erties 



Keywords 

Fokker-Planck equation, continuum mechanics, incompressible fluid 



1 Introduction 

The object of our considerations is a special case of Fokker - Planck equation, which describes evolution 
of 3D continuum of non-interacting particles imbedded in a dense medium without outer forces. The in- 
teraction between particles and medium combines diffusion in physical space and velocities space. Classic 
Fokker - Planck equation contains only two forces, which act on particles: damping force proportional to 
velocity and random force. 

The aim of this work is to derive equation for the case of constrained movement - the space density of 
particles remains constant. As a result, our equation contains the third force - constraint reaction. This 
is equivalent to hydrostatic pressure. 

Our starting point is the classical Fokker - Planck equation: 




(1) 



where 



n = n(t,xi,X2,X3,t;i,t;2,W3) - density; 
t - time variable; 
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Xi,X2,xs - space coordinates; 

vi,V2,vs - velocities; 

a - coefficient of damping; 

k - coefficient of diffusion. 

In the following section we add to this equation terms due to reaction of constraints and thus de- 
rive the Fokker - Planck equation for incompressible fluid. Then we rewrite the equation for curvilinear 
coordinates. In the fourth section we obtain a stationary solution of Fokker - Planck equation for in- 
compressible fluid and in the fifth section we derive linearized equation. In the last section we introduce 
symmetries of Fokker - Planck equation for incompressible fluid. 

2 Derivation of Fokker - Planck equation for incompressible fluid 

In our previous work (see ref. [1-2]) we deducted equation (1) from the kinematic equation (see below). 
The crucial point is expression for the force, which acts upon particles. We represented the force as a 
sum of random force Ri and friction force, which is proportional to the velocity: 

Qi = Ri- avi. (2) 
We calculated random force i?j from the law 

„ , dn 

nR, = -k—^. (3) 

(2-3) give exphcit expressions for the forces. As we show in [2], (2-3) imply equation (1). 

Now we try to get expression for reaction of constraint. We need some additional principle for this 
purpose. In the following we accept the Gauss principle of least constraint (see [3]). 

Gauss principle reads for discrete system of particles: 




where 

Xr,yr,Zr - material point coordinates ; 
Xr,yr,Zr - accelerations ; 

Xr, Yr, Zr - Components of outer forces, acting on this point. They are dependent on coordinates Xr, Vr, Zr 
and velocities Xr,yriZr- 

Besides X^, YJ,, reactions of constraints are also to be considered, though they are not given explic- 
itly. They appear as Lagrange multipliers of constrains. 

Let us suppose, that coordinates Xr,yr,Zr and velocities kinematically possible, that 

is they satisfy constraints (holonomic or nonholonomic) . Then effective accelerations Xr,yr,Zr provide 
minimum of G, comparing with all kinematically possible accelerations. 
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The derivation of reactions of constrains for the simple case of discrete system of particles is straight- 
forward (see [3]). We deal with more complicated case of continuum with dispersed velocities - many 
particles with different velocities are simultaneously placed in each point of space. 

But first of all we examine first intermediate case of continuum with uniquely defined velocities field, 
which is simpler. For this case the Gauss principle was used in [4], [5]. 

For such a continuum (fluid) we have the following expression for accelerations: 

dui dui 

In this way we get the following expression for Gauss functional: 

G = [ ^ (am - -Fm^ (am - -Fm] dV. (6) 
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Condition of incompressibility of the velocities field Ui is 



0. (7) 



dxk 

Differentiating (7) on time, we get expression for constraint on accelerations field: 



dxkdt 

This implies the following constraint on accelerations a*: 



^'"^ =0. (8) 



dak diik , Of dum\ duk dum 



dxk dxk dxm \ dxk J dxm dxk 

To take in account (9) we use Lagrange multipliers method. This gives following modified expression 
for Gauss functional: 



G 



2 V P J \ P J \OXk OXm OXk 



dV. (10) 



where p is Lagrange multiplier for constraint (9) . We seek the extreme of functional (10) with Uj as 
unknown variables. 

Variational Euler equations for functional (10) are 

Um + n^--Fm)+^=0. (11) 
dXi p J dXm 

They are identical to standard Euler equations for incompressible fluid. We see, that Lagrange's 
multiplier p is pressure and 1/p {dp/dx^) is expression for reaction of incompressibility constraint. 

We conclude, that as for get Fokker - Planck equation for incompressible fluid we should: 

- to add to system additional equation - expression for incompressibility constraint; 
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- to add to all forces, acting on particle (see (2)), gradient of the pressure as reaction of ideal constraint 



Q, = R^-^^-avi. (12) 

In this way we get the following system with one integral constraint expression and one differential 
equation. 



p = n dvidv2dvs = const, (13) 
Jv 



where p is constant density. 



dn ^ dn d{vjn) 1 dn dp ^ d'^n 
dt dxk dvj pdvjdxj dvjdvj 



(13) imply equations {dp/dxi) = 0. 

This derivation method is insufficient, because we obtained expression (12) from the wrong model. 
Now we return to the case of continuum with dispersed velocities - there exist many particles with different 

velocities in each point of space. 

Kinematic of single particle is described by differential equations: 

(la) 

where 

Xk - space coordinates; 
Vk - velocities of particles; 
bk - accelerations of particles. 

Evolution of particles density n satisfies the following equation: 

dn ^ djnvk) ^ djnbk) ^ ^ 
dt dxk dvk 

If accelerations are independent from n, (15) and (16) are equations of characteristics of (17). 
The dynamics of the system is equation of movement of particles - the second Newton's law: 

b, = 1q,; (18) 
where m - particles mass, the force Qj is the same as in (2-3), (12). 
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We accept (13) as constraint on n and try to get expression for constraints on accelerations bj. (13) 
implies 

/dn 
^ —dvidv2dv3 = 0, (19) 

f dn 

/ —dvidv2dv3 = 0, (20) 
Jv oxj 

This means, that averages of n derivatives on space and time coordinates are zero. 
Let us integrate (17) on velocities and use (20). We get identity: 

f ^ipildv,dv2dvs = 0. (21) 
Jv oxk 

But we know, that average velocity is equal to 

puk = / nvkdvidv2dv3. (22) 
Jv 

Therefore (22) is similar to (7) - the divergence of average velocity field is zero. We can write this in 
the following way: 

f dn 

/ ^ — Vkdvidv2dvs = 0. (23) 
Jv oxk 

This identity concerns the first moments of derivatives of n on space coordinates. They are not zero, 
but their sum is zero. 

Let us multiply (17) by vj and integrate the result by parts: 



djnvj) ^ djnvkVj) ^ djnbkVj) __ 
dt dxk dvk 



Let us integrate (25) on velocities. We get expression for average accelerations (compare this with 

(5)) 

p{bj)avg = J^nbjdvidv2dv3 = ( ^^^J'^'* + ^^^Xfc^^O *^25) 



Let us denote 



Jkj = / nvkVjdvidv2dv3 = pukUj - Ukj, (26) 

J 



where cTfcj - components of stresses tensor. 
So (26) reads 



P[bj)avg- + g dXk~^' dXk- ^ ^ 



5 



This is very remarkable, that average acceleration is not equal to acceleration of average motion and 
contain additional term - divergence of stresses tensor. This follows from the fact, that accelerations are 
quadratic on velocities. 

We shall use (27) to get constraint expression. 

Integrate Newton's law (18) and get: 



Pibj)avg = (Qj) 

(27) and (28) imply equations of movement : 



avg 



(28) 



dxk 



(29) 



The Gauss functional without constraints has the following form 



G = dxidx2dx3 
Jx 



Jv ^ \ rn 



h- —Qk] dvidv2dv3. 
m 



(30) 



To get expression for constraint on accelerations we take divergence of (27) (compare with (9)). 

d , , , , , , f d f d{nvj) d{nvkVj) 



I —— (nbj) dvidv2dv3 = [ -— 
Jv oxj Jv oxj 



dt 



+ 



dxi 



dvidv2dv3. 



(31) 



The first term in the RHS of (31) is zero according to (21). The rest terms give us the desired 
constraint expression: 
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n 



dxjdxk 



VkVj 



dvidv2dv3 = 0. 



(32) 



We multiply (32) by p/p and add the result to (31). Gauss functional with constraint reads now 
(compare with (10)): 



G = [ dxidx2dx3 [ (77 (bfc - —Q 
Jx 1 2 V m 



k j \bk Qk 

m 



(33) 



p{xi 



_d_ 

dx-. 



(nbj) 



n 



dxjdxk 



VkVj 



dvidv2dv3. 



We seek the extreme of functional (33) with (nbj) as unknown variables. The Euler's equation for 
(33) is: 



n I dp 

nbk = —Qk n-—. 

m p oxk 



(34) 



This means, that reaction of constrains is really equivalent to pressure field and we return to (12) and 
(14) once again. 
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3 Fokker - Planck equation for incompressible fluid in curvilinear co- 
ordinates 



Now let us consider the form of (13-14) system in the curvihnear coordinates. We need not perform all 
calculations here, because we can refer for details to our previous works [1-2]. Therefore we give here 
only results of calculations. 

Let us denote gmn - the covariant components of metric tensor and g"^"' - the contravariant components 
of metric tensor, g = det\gij\ = l/det\g''^ . 

Christoffel's symbol is called 

J. _ 1 / dgnp dgpm dg 

mn \ /nr\ 

J- p,mn — o I 'aZ:^ + "aZ^ KZT I ■ l^"^-* 



2 V^X"^ dx^ dxP . 
and 

^mn ~ • (36) 

Covariant components of velocity vector are Vi, contravariant components of velocity vector arc v^. 

Using these definitions, we can write equations (13-14) in curvilinear coordinates with contravariant 
velocities as independent variables, as 

J dv^dv'^dv^ = 1. (37) 

^ + -'^^ - - « -'^^ - 3 « - (38) 

„rnn 



For details we refer to our works [1-2] . 

For example, the system (37-38) in spherical coordinates reads: 

/ n dv^dv^dv^ = „ ^ (39) 
Jv sm{d) 

dn -,dn r,dn ndn 

+r Uv'' + siu\eyv^) 1^ + f sin(^) cosiO^v^ - -v^A |^ - 2 f + ^^^v\A 
^ ^ ov^ \ r J ov^ \ sin[6) J ov^ 

f I dn 2 3 \ 2 1 ( dn dp 1 dn dp 1 dn dp\ 

\ dv^ dv'^ dv^ J pKdv^dr r'^ dv"^ dO sin^ {6) dv^ d4> J 

^ / d'^n 1 d'^n 1 d'^n \ 

\dv^dv^ r"^ dv'^dv'^ r'^ sm^(9) dv^dv^ j 

The system (13-14) in curvilinear coordinates and with covariant velocities as independent variables, 

is 
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/ n — dvidv2dv3 = p. (41) 
Jv ^Jg 

dn , mk 9n i dn dn 

-Qi + a ^-d^k+^kifv,v,^^-av,g^^-San- (42) 

1 dn dp , d^n 
- k 9ik- 



p OVk ox" OVlOVk 

For example, the system (41-42) in spherical coordinates has the form: 

/ n dvidv2dvs = r^sm(9)p. (43) 
Jv 

dn dn V2 dn dn , , 

'm'^'^^d^^^de^ r2sin2(0) d^^ ^ ' 

+J_ (v V + ^^^^ \ — + '^°^(^)'"3^3 dn 
y sin^(6') J dvi r2sin^(^) dv2 

f dn dn dn\ 1 f dn dp dn dp dn dp 
—a vi- h ^27^ 1- '"37^ — —San tttw + "TT"^ + 



dvi dv2 dv^J p\dv^dr dv"^ dO dv^ d(j) 

, / (9^n 9 d'^n o . 9,^, d'^n \ 



dvidvi dv2dv2 dvsdvs I 

For orthogonal coordinates the diagonal components of metric tensor are expressed as squares of Hi 
- Lame coefficients. All off-diagonal components are zero. The Christoffel's symbol components can be 
expressed as derivatives of the Lame coefficients. 

The system (13-14) in curvilinear coordinates and with physical velocities as independent variables, 

is 

/ n dw^dw^dw^ = p. (45) 
Jv 

dn dn dn w'^ / g'^^s kd^k\ fc ^ /^g\ 

dt Hk dx^ dw^ HsHk \ dx^ dx^ ) dw^ 

11 dn dp _ ^ d^n / j _ ^\ 
p Hk dw^ dx^ dw^dw^ ' 

For example, the system (45-46) in spherical coordinates has the form: 

/ n dw^dvo^dw^ = p. (47) 
Jv 

dn 1 dn w"^ dn w'^ dn , , 

— + — — + ^-— — + 48 
dt dr r dO r sm(^) dcp 

I/99 qo\9n If cos(6) 00 1 n\ dn 1 / 1 q cos(0) 9 q\ dn 

H — iw'^w'^ + W^W^ -— j- + - ^W^W^ - w^w^ I ( w^W^ -\ j^w w ' 

r \ ) dw^ r \i 



dw^ r \sin(^) / dw'^ r \ sm{9) J dw^ 
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f -I dn dn n dn \ 
- a \ w ——r + w -— ^ + w - 3 a n- 

V ow'- ow^ ow^ / 

1 f dn dp \ dn dp 1 dn dp\ ^ ( d'^n d'^n d'^n \ 

p\dv^dr r dv'^ do r sm{9) dv^ d(p J \dw^dw^ dw^dw^ dw^dw^ j 



4 Stationary solution of Fokker - Planck equation for incompressible 
fluid 

Only in this section we discuss the case of nonzero force field. In this case equation (14) reads 

dn dn d . , 1 9n / dp\ , d'^n , ^ , 

where Fi{x) are components of force acting on particle. For non potential forces the stationary solution 
of Fokker - Planck equation does not exist. For potential forces we have the following expressions for 
force components: 

and (14) reads 

dn ^ dn d . ^ ^ (p- \ ^ C51'\ 

dt ■' dxj dvj ■' dvj dxj dvjdvj ' 

where Pi{x) - potential function. 

Stationary solution of usual Fokker - Planck equation has the form n = m{x) s{v). For the case of 
incompressible fluid must be m{x) = 1 because of condition (13). Therefore n = n{v). 

To kill all the terms with cross products of derivatives on space coordinates and velocities in (51), we 
substitute: 

p = const — Pi. (52) 

This is the Pascal's law. 

The rest of terms in (51), depending only on velocities, is the divergence of some current in velocities 
space. For the true static solution all components of the current must be zero: 

dn 

—avjn = k— — . (53) 
dvj 

So n has Maxwell distribution : 

(54) 



a 

V.\V.\ 

2k 



n = no exp 

The value of no constant factor we find from condition (13) 
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no = p 



a 



y/2 



(55) 



a 



(56) 



5 Linearization of Fokker - Planck equation for incompressible fluid 

We seek solution of (13-14) as a sum of stationary solution and perturbation term: 



p = ep. 



(57) 



a \3/2 



a 



+ en. 



(Pi = here). 

Substitute these expressions to (13-14) and drop terms of more than first degree on e. 



(58) 




The system (59-60) though linear, is not trivial. We shall discuss it's solution later. 

6 Symmetries of Fokker - Planck equation for incompressible fluid 

We presented symmetries of standard Fokker - Planck equation (1) in our previous work [6] with all 
calculations details. In this section we present only result of calculations for Fokker - Planck equation for 
incompressible fluid. 

Two equations of the (13-14) have different nature: (13) is integral equation and (14) is differential. To 
obtain symmetries of this system we perform two-step process: 1) we determine symmetries of differential 
equation (14) and 2) we calculate actions of these symmetries operators upon (13) and drop all invalid 

operators, for which action does not vanish. 

Wc omit tedious details of step 1) and give below only resulting expressions. Expressions for variations 

of variables are: 



Sx = r2Z - ray + 3 Cie°^/2*a; + /i(t); 



(61) 
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5y = nx - + 3 Cie°*/2y ^ ^^^^y^ 
5z = ny -r2X + 3 Cie"*/^^ + fsit). 
where r\,r2,rz-,C\ - arbitrary constant coefficients; /i,/2,/3 - arbitrary functions of argument t. 

5u = Cie^^/^x + de^^/^u + r2W - + f[{t)- 
5v = Cie'^^/^y + Cie"*/^!; + r^u - nw + f^{t); 
Sw = Cie"*/2z + Cie°*/2i(; + nv - r2U + /^(t); 



(62) 



St = -Cie"*/2 + C2. 
a 



(63) 



Sn = /4(ne3°*) e'^''* - VlCxe"-^l''n. 



(64) 



5v = {ccf[ - f1)x + {af'2 - f!i)y + {af^ - f'i)z+ 
+la^C,e'^'l\x^ + / + z^) + 2Cie"*/2p + f,{t). 



(65) 



where C2 - one more arbitrary constant coefficient; /4, /s - two more arbitrary functions of argument t. 

These expressions lead to the following expressions for symmetries operators: 
- operator (rather exotic) associated with Ci: 



'^d /3 9, 9 9 9x ^ \ ^ 45 



/3 \9/3 \9/3 \ d d d d 

time shift operator associated with C2: 



(66) 



■"2 



at' 



- three operators of shifts along x,y,z axes associated with /i,/2,/3 functions: 



(67) 



^4 = («/^-i^0.| + /^w| + /.(4; 
-a = M-/^>| + /^W^ + /3(4; 



(68) 

(69) 
(70) 
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- one more exotic operator associated with /4: 



3at\ „— 3at 



V6 = f4{ne'n e 



dn 



(71) 



- operator, which states, that we can freely choose arbitrary additive pressure at each moment of time, 
associated with /s: 



- three time independent rotations: 



■U8 = 


d 

dw 


d 
dv 


d 


d 


Vg = 


d 
du 


d 

dw 


d 
dx 


d 

dz' 


VlO = 


d 
dv 


d 

du 


d 
dy 


d 

-y^x 



Now expression for variation of (13) is: 

5n + n 



V 



d6u d6v d6w \ 

+ ^ + -r— I 



du dv dw J 



dvidv2dv3 = 0, 



Substitute (66-75) to (76) and get: 
Jv L 



dvidv2dv^ = 0, 



This means 



and 



Ci = 0; 



/4 = 0. 



(72) 



(73) 
(74) 
(75) 



(76) 



(77) 



(78) 



(79) 



Therefore symmetries of Fokker - Planck equation for incompressible fluid build the subgroup of 
(66-75) with excluded exotic symmetries vi and vq. 

This list of symmetries is the same as the Navier - Stokes equations symmetries list (see ref. [7]), 
with exception of scaling symmetries. Our system has no scaling symmetries, because it contain additive 
frictional term. 
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7 Symmetries of linearized equations 



Similar to previous section, we present only result of symmetries calculations without calculations details. 

We use the same two-step process. As a result of the first step we get following expressions for 
variations of variables: 

5n = Cm + n, (80) 

6p = Cip + p, (81) 
where Ci - arbitrary constant coefficient, n, p - arbitrary solutions of linearized system (59-60) itself; 

6t = C2, (82) 
where C2 - one more arbitrary constant coefficient; 

Sx = r2Z - rsy + C3, (83) 
Sy = r2,x -riz + Ci, 
Sz = riy-r2X + C5, 

where ri, r2, ra, C3, C4, C5 - arbitrary constant coefficients; 

Su = r2W — r^v, (84) 

5v = r^u — riw, 
6w = riv — r2U. 

These expressions lead to the following expressions for symmetries operators: 

- scaling operator associated with Ci: 

"^ = "l^+^|' ^^^^ 

- time shift operator associated with C2: 

V. = |; (86) 

- three operators of shifts along x, y, z axes associated with C3, C4, C5 : 

»3 = I:; (87) 
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- three rotations: 





d 




u 

dw 


u u 
-w— + y— 
ov oz 


f) 
(J 

dy 


d 

du 


d d 
ow ox 


d 


d 

dv 


du dy 


d 



- infinite subgroup of finear system (59-60) solutions: 



_d _d 
vq = n— +p^. 
dn dp 



Now we go to second step. Expression for variation of (59) is: 

dvidv2dv3 = 0, 



V 



( d5u d8v dSw 
on + n — 1- — h 



V du dv dw 



Substitute (80-84) to (94) and get: 



/ (Cin) dvidv2dvs = 0, 
Jv 



This means 



Ci = 0. 



(89) 

(90) 
(91) 
(92) 



(93) 



(94) 



(95) 



(96) 



Therefore symmetries of hnearized Fokker - Planck equation for incompressible fluid build the sub- 
group of (85-93) with excluded scaling symmetry vi. 

Let us give some small examples of invariant solutions. 

Rotations (90-92) have following invariants: 



r = ^x'^ + y"^ + z'^; 



U = Vu"^ + v'^ + w'^; 
Ur = xu + yv + zw. 



(97) 

(98) 
(99) 



According to PDE symmetries theory (ref. [7]), we should seek solutions of (59-60) of the following 
form: 
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n = n{t, r, U, 11^); 



(100) 
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DISCUSSION 



In this article we deducted the Fokker - Planck equation for incompressible fluid from the Gauss principle. 
This method can be useful for investigation of another cases of constrained movement. To promote this 
task we give three slightly different forms of equations for use of curvilinear coordinates. 

Unfortunately, the actual solution of equations is a complicated task, because we move from the model 
with totally independent particles to model with interaction. This model is by necessity nonlinear. To 
make first steps to solution we use perturbation method. We construct linearized equations, but do not 
try to solve them at the moment - this is a task for another work. 

It is of some interest to obtain symmetries group for our system. We find the group and compare it with 
symmetries of Navier - Stokes equations. One symmetry is missing in our case as a result of friction force 
presence. 

We derive symmetries of linearized equations also and give some examples of invariant solutions. We see, 
that symmetries group of linearized system is not rich enough and is hardly usable to obtain physically 
interesting solutions. 
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